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Abstract 

p 

The India-based Neutrino Observatory (INO) will host a 50 kt magnetized iron calorime- 
ter (ICAL@INO) for the study of atmospheric neutrinos. Using the detector resolutions and 
efficiencies obtained by the INO collaboration from a full-detector GEANT4-based simula- 
I tion, we determine the reach of this experiment for the measurement of the atmospheric 

neutrino mixing parameters (sin 2 623 and |Am| 2 |)- We also explore the sensitivity of this 
experiment to the deviation of 823 from maximal mixing, and its octant. 
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Neutrino flavor oscillations have been established beyond any doubt by a series of outstanding 
results from atmospheric [I], solar j2j, reactor [3jj6] and accelerator [7jj9] neutrino experiments. 
Neutrino flavor oscillations require neutrinos to be massive and mixed, and they have thus pro- 
vided the first unambiguous hint for physics beyond the standard model of elementary particles. 
The neutrino mixing matrix, called C/pmns (l0, 11 , can be parameterized in terms of three mix- 



5_i 

ing angles #12, #13, #23) arid a charge-parity violating (CP) phase 5cp- In addition, if neutrinos 
were Majorana particles, we would have two Majorana CP phases a\ and 02 as well. The 
frequencies of neutrino oscillations are governed by two mass squared differences, Am^ and 
Am^, where we define Am 2 ^ = m 2 — m 2 . The so-called solar neutrino oscillation parameters 
612 and Am^ have been measured from the combined analysis of the KamLAND reactor v e 
data and the solar neutrino data. The so-called atmospheric neutrino oscillation parameters #23 
and I Am^ I are mostly constrained by the Super-Kamiokande (SK) atmospheric and MINOS f« 
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disappearance data. The third and the last mixing angle #13 is the latest to be measured by a 
series of accelerator and reactor experiments. After decades of speculation on whether #13 was 
zero, data from these experiments have revealed that the value of $13 is not only non-zero, it is 
in fact just below the previous upper bound [12] from the Chooz experiment. Accelerator-based 
neutrino experiments T2K and MINOS have both observed v e appearance events from a beam 
of Vn that indicates a non-zero value of #13. The short baseline reactor neutrino experiments 
Daya Bay, RENO and Double Chooz have excluded #13 = at 5.2a, 4.9a and 3.1a respectively 
from disappearance. Their best fit values are sin 2 2^13 = 0.092 ± 0.016(stat) ± 0.005(syst) [41, 
sin 2 26 13 = 0.113±0.013(stat)±0.019(syst) [5] and sin 2 2% = 0.109±0.030(stat)±0.025(syst) [6], 
respectively. We summarize our current understanding of the neutrino oscillation parameters in 
TableE 



Parameter 


Best Fit Value 


3a Ranges 


sin 2 #12 


0.307 


0.259-0.359 


sin 2 6>23 


0.386 


0.331-0.637 (NH) 




0.392 


0.335-0.663 (IH) 


sin 2 6*i3 


0.0241 


0.0169-0.0313 (NH) 




0.0244 


0.0171-0.0315 (IH) 


Am 2 ,! (eV 2 ) 


7.54 x 10~ 5 


6.99-8.18 x 10~ 5 


lAmfil (eV 2 ) 


2.43 x 10~ 3 


2.19-2.62 x 10~ 3 (NH) 




2.42 x 10~ 3 


2.17-2.61 x 10" 3 (IH) 



Table 1: A summary of the current values of the neutrino oscillation parameters. The values 
are taken from |13| . NH and IH indicate the limits applicable if the mass hierarchy is normal 
and inverted, respectively. 

While the above experiments continue to improve the precision on the mixing angle #13, the 
focus has now shifted to the determination of the other unknown parameters in the neutrino 
sector. These include the CP violating phase 5cp and the sign of Am^. Though the current 



global analyses 13, 14 have started to provide some hints about the value of Sep, these are still 
early days, and better data from dedicated experiments would probably be required before one 
could make any definitive statement on whether CP is violated in the lepton sector as well. The 
answer to this question would have far-reaching implications, as a positive answer would lend 
support to the idea of baryogenesis via leptogenesis during the early universe. 

Determining the sign of Am]j is also of crucial importance, since its knowledge is essential 
for constructing the mass spectrum of neutrinos. This sign can be extracted via the detection of 
earth matter effects in atmospheric and accelerator-based neutrino beams. Proposed detectors 
aiming to observe earth matter effects in atmospheric neutrinos include the Iron CALorimeter 
(ICAL) at the India-based Neutrino Observatory (INO) [15], megaton water Cerenkov detectors 



such as the Hyper-Kamiokande (HK) [16], large liquid argon detectors 17 [18] and a gigaton- 
class ice detector such as the Precision IceCube Next Generation Upgrade (PINGU) [19 20 . 
The accelerator beam neutrino experiment NOz/A [21J, which will start operating soon, will also 
be sensitive to the mass hierarchy. 

The other neutrino parameter that is yet to be measured is the absolute neutrino mass scale, 



on which currently we have only upper bounds from cosmological data |22-24 , from neutrinoless 
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double beta decay 25,26 , and from tritium beta decay experiments 27,28 . Finally, the neutri- 



noless double beta decay experiments will also answer the most fundamental question - whether 
neutrinos are Dirac or Majorana particles. If neutrinos are indeed Majorana, the neutrinoless 
double beta decay experiments might also have a chance to shed light on the Majorana phases 
in the far future. 

INO is an experiment |15| proposed to study atmospheric neutrinos, employing the ICAL 
that can distinguish between /x + and fx" , and will have good energy as well as direction resolution 
for muons. The magnetic field will enable it to distinguish between neutrinos and anti-neutrinos, 
hence allowing identification of the neutrino mass hierarchy. It will be located in the Bodi West 
Hills in the Theni district of Tamil Nadu in Southern India. The ICAL cavern will be located 
under a 1589 m high mountain peak. This peak provides a minimum rock cover of about 1 km in 
all directions to reduce the cosmic muon background. The detector is designed to consist of 150 
alternate layers of 5.6 cm thick iron plates and glass Resistive Plate Chambers (RPCs) stacked 
on top of each other. The total mass of the detector will be about 50 kt. The iron plates act 
as the target mass for neutrino interactions and the RPCs are the active detector elements that 
give the trajectory of the muons and hadrons passing through them. Iron will be magnetized to 
a field of about 1.3-1.5 Tesla. 

The main physics objective of ICAL is the determination of the neutrino mass hierarchy 
through the observation of earth matter effects in atmospheric neutrinos. In a recent anal- 
ysis 29 performed by members of the INO collaboration, it was shown that data from this 



experiment along with that from the reactor and long baseline experiments T2K and NOi^A 
could determine the neutrino mass hierarchy at the 2.2cr — 5.5cr C.L., depending on the true 
values of the parameters #13, #23 and <5cP) with 50 kt x 10 years of exposure. 

In this paper, we explore in detail the potential for measuring the neutrino parameters #23 arid 
|Am| 2 | in the ICAL@INO experiment using atmospheric neutrinos. The precision on both these 
parameters is expected to improve from data coming from the currently operating and soon-to- 
start experiments using accelerator-based neutrino beams (MINOS, T2K and NOz^A) as well as 
atmospheric neutrinos (at Super-Kamiokande and IceCube Deep Core). Here we compare the 
potential of the ICAL@INO atmospheric neutrino experiment vis-a-vis the expected sensitivity 
of these complementary experiments. For our simulation of atmospheric neutrino events in 
ICAL@INO we use the software tailored and developed by the INO collaboration and present 
results obtained using the detector resolutions and efficiencies obtained from detailed simulations 



of the ICAL detector 30,31 . We restrict ourselves to using the information about the energy 
and direction of the muons produced in the detector. The results reported here are a part of 
the simulation and physics analysis being performed by the INO collaboration towards the full 
understanding of the capabilities of this experiment and a realistic study of its sensitivity reach. 
(Studies that explore the reach of a magnetized iron calorimeter for neutrino mixing have been 



performed in the past 32 33 , however this is the first time that we use the realistic ICAL 
detector simulations in the analysis.) 

The paper is organized as follows. In section 2 we describe the detector simulation proce- 
dure for the ICAL detector, and the detector response for muons obtained from the GEANT4 
simulations which are used in this work. In section 3, we outline our analysis procedure. In 
section 4, we discuss the results of our analysis, presenting the reach of the ICAL@INO for #23 
and IA777.32I. Section 5 summarizes our findings. 



3 



2 Simulation Framework and the Detector Response for Muons 



Neutrinos interact in the ICAL detector with nucleons and electrons, giving rise to leptons and 
hadrons. As the outgoing lepton and hadrons propagate through the detector volume, they 
cross one or more RPC layers to give hit points in the RPCs. For every hit point, the (x, y) 
position, time of hit and the RPC layer number are recorded. The information contained in 
the hit points is used to reconstruct the energy and direction of the particles. A muon passing 
through the ICAL detector gives hit points in several layers, which can be cleanly joined to 
form a track. In contrast, a pion or other hadron passing though the detector typically gives 
rise to a hadron shower with multiple hit points in a single layer which cannot be joined to 
form a track. The ICAL detector is optimized primarily to measure the muon momentum with 
a good precision. The muon charge identification efficiency is also high. The hadron energy 
and direction can be calibrated based on the number of hit points, and by combining the muon 
and hadron information, the neutrino momentum can be reconstructed, albeit with a coarser 
precision. 



The detector simulation procedure is as follows. A GEANT4 134,35 -based detector simula- 



tion package for ICAL has been developed. The secondary particles are first propogated inside 
the detector using the GEANT4 code, taking into account the energy losses, multiple scatter- 
ing, and magnetic field. This gives the hit positions of the particles in the active detectors, i.e. 
RPCs. Next, this output is digitized and used as an input for the reconstruction of the particle 
energy and direction. 

The reconstruction of the momentum of a particle is done in two steps. First, the topology of 
the digitized hit points is analyzed to find a possible muon track, and/or to find a hadron shower. 
Next, if the hit points are found to form a muon track, a Kalman Filter-based algorithm is used 
to reconstruct the initial energy, direction, and the vertex position of the muon. In case the hit 
points are identified as a hadron shower, the energy of the hadron is estimated by calibrating it 
against the number of hit points. 

The muon response of the ICAL detector is parametrized in terms of the energy resolution 
(<xe), direction resolution (a cos g), reconstruction efficiency (e#) and the charge identification 
efficiency (ec)- All these quantities are estimated as functions of the true muon energy (E^) and 
direction (cos 9^), for and fi + . Here 8^ is the zenith angle for the muon. Figure [l] shows the 
muon energy and direction resolutions, as well as the reconstruction and charge identification 
(CID) efficiencies for for some specific muon zenith angles. The numbers for fi + are virtually 
identical. We note here that the reconstructed energy distribution of the muon is fitted with 
the Normal distribution function for > 1 GeV, whereas for E^ < 1 GeV it is fitted to the 
Landau distribution function. In the present work, we have carried out the oscillation analysis 
with muon energy and direction (E^, cosO^) as observables using simulated data. 

The analysis in this paper is based only on the energy, direction, and charge of the muon, 
for the charged-current (CC) events. The inclusion of information on hadron energy and 
direction is beyond the scope of this paper. The detector response to muons used in this work 
comes from the first set of simulations done with the ICAL code. These simulations are ongoing 
and are being fine-tuned. Hence the values of the resolution functions and efficiencies are likely 
to evolve along with the simulations. 
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Figure 1: The energy resolution (a), cosO^ resolution (b), the reconstruction efficiency (c), and 
the charge identification efficiencies (d), for \x~ as a function of the true muon energy and for 
three cases of true muon zenith angle. The red, green and blue lines are for the zenith angle 



bins with cos 6„ 



-0.35 ± 0.05, -0.55 ± 0.05, -0.85 ± 0.05, respectively. 



3 Oscillation Analysis Procedure 

The physics analysis of atmospheric neutrino events requires simulations which can be broadly 
classified into four steps : (i) neutrino event generation, (ii) inclusion of the oscillation effects, 
(iii) folding in the detector response, and (iv) the \ 2 analysis. These steps are described in detail 
in the following subsections. 

3.1 Event Generation 

We use the neutrino event generator NUANCE (version 3.5) [36] to generate the neutrino inter- 
actions. The atmospheric neutrino fluxes provided by Honda et al. |37| at the Super Kamiokande 
location are usedj^] The ICAL detector composition and geometry are specified as an input. NU- 
ANCE incorporates the differential cross section for charged-current (CC) and neutral-current 

"The neutrino fluxes at Theni have not yet been incorporated in the NUANCE event generator. 
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(NC) interactions for all nuclear constituents of the materials used in the detector, for all neu- 
trino flavors, and for all possible interaction processes. The quasi-elastic (QE), resonance (RS), 
deep inelastic (DIS), coherent (CO) and diffractive (DF) scattering processes are all included. 
The neutrino fluxes are then multiplied with the interaction cross-sections to calculate the event 
rates for all scattering processes and possible target nuclei. Event kinematics are generated 
based on the differential cross- sections. 

The NUANCE output consists of the 4-momentum (p^) of the initial, intermediate and the 
final state particles for each event. To reduce the Monte Carlo (MC) fluctuations in the number 
of events given by NUANCE, we generate a very large number of neutrino interactions (an 
exposure of 50 kt x 1000 years) and scale it down to the desired exposure for the \ 2 analysis. 



3.2 Inclusion of Oscillations 



The total number of Un events coming from the 



Vn and the v e 



channels is given as 



d 2 N \ d 2 $ Vu <f$ v 

—— —r = N T xN D x a v x p ___Ji__ + p « 

dE v d[cosB u ) ^ dE v a(cos U v ) dE u d(cosO, 



(1) 



where Nt is the exposure time and Njj is the number of targets in the detector. Here <& v and 
<& Ue are the fluxes of and u e respectively, and P a $ is the u a — > up oscillation probability. 
Since an exposure of 1000 years has to be taken to reduce the MC fluctuations, generating 
events using NUANCE for each set of oscillation parameters is extremely time consuming and 
practically impossible within a reasonable time frame. Therefore, we generate events for 1000 
years of exposure using NUANCE only once and employ the event re- weighting method described 
below to take care of neutrino oscillations on this event sample for any given set of oscillation 
parameters. 

To get the number of \x~ events from z/„ that have survived oscillations, we generate the 
Vn interactions from NUANCE using the un-oscillated flux in the neutrino energy range 
0.5 — 100 GeV and all neutrino zenith angles 9 U . For a given cos 6 U , the path travelled between 
the production point and the detector is 



L = \J (R + L ) 2 - (R sin9 u ) 2 - R cos9 h 



(2) 



where R is the radius of the earth (6378 Km) and Lq is the average height of the atmospheric 
neutrino production, taken here to be 15 km. The survival probability P w for the neutrino 
energy E v and zenith angle 9 U for a given event is calculated using the oscillation probability 
code included in NUANCE [36]. The details of this oscillation probability calculation are as 
described in [?]. We next impose the event re- weighting algorithm as follows. To decide whether 
an un-oscillated event survives the oscillations to be detected uniform random 

number r is generated between and 1. If P^ e < r < P^ e + Pnui we keep this event as a z/« 
event. Otherwise this v„ is considered to have oscillated into a different flavor. 

Similarly, to get the \x~ events from atmospheric v e converted to due to flavor oscillations, 
we generate the Vn interactions from NUANCE using the un-oscillated v e flux in the neutrino 
energy range 0.5 — 100 GeV and all neutrino zenith angles U . This then has to be folded in with 
the conversion probability P e „. For that we again use the event re- weighting algorithm where 
for each event a random number r is generated between and 1. If r < P efl we label the event 
as a event, otherwise we discard it. 
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The Un events from Vn —> Vn as well as v e — > oscillation channels are added to form the 
[x~ event sample. Similarly, we form the y + event sample as well. Next, we read the muon 
momentum of all the events that we have collected after oscillations and bin them according to 
energy (En) and direction (cos0^) of the muons. We bin the data in the energy range En = 0.5 
- 15.5 GeV (300 bins) and cos0 M in the range —1 to +1 (20 bins) because our analysis threshold 
is 0.8 GeV and the muon reconstruction efficiency is very low below 0.5 GeV. We keep track of 
the numbers of yT and /x + events separately. At this stage, we have the distribution of muon 
events in terms of their "true" energy (En) and the cosine of the zenith angle (cos 6^). 

Figure [2] shows the zenith angle distribution of yT events in the energy bin 2 — 3 GeV before 
and after invoking oscillations using the re-weighting algorithm described above. We use the 
oscillation parameters described in Tableland take the exposure to be 50 kt x 10 years. 



Parameter 


sin 2 2#i2 


sin 2 20 23 


sin 2 20i3 


Am 2 -! (eV 2 ) 


|Am 2 2 | (eV 2 ) 


^CP 


Hierarchy 


True Value 


0.86 


1.0 


0.113 


7.6 x 10" 5 


2.242 x 10" a 


0.0 


Normal 



Table 2: True values of the neutrino oscillation parameters used in this paper 
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— if— With oscillations ~ 

2 GeV < E„ < 3GeV I _Z 

50 kt x 1 years ~ 
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-1 -0.5 0.5 1 

cos 9^ 

Figure 2: Zenith angle distribution of fi~ events for the bin 2 GeV < En < 3 GeV, without and 
with flavor oscillations. The detector efficiencies have not been included here. The error bars 
shown here are statistical. 

3.3 Folding in the Detector Response 

Next, we fold in the detector response [30] to obtain the measured distribution of muons. We ap- 
ply the reconstruction efficiency (e_R- ) for yT by multiplying the number of events in a given true 
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Figure 3: Zenith angle distribution of oscillated \i events for the bin 2 GeV < En < 3 GeV, 
after taking into account detector efficiencies. The error bars shown here are statistical. 



bin with the corresponding reconstruction efficiency: 

x N* ue ' 



energy (E^) and true zenith angle (cos 6^ 

iV cos 0„) = e R - (E^, cos x N^{E^ cos 9 U ) , (3) 

where iV*™ e is the number of events in a given (Ep, cos#^) bin. Exactly the same procedure 
is used for determining the /i + events. The CID efficiency (ec- for [i~ and ec+ for n + event 
sample) is next applied as follows: 



+ (1 - e c +) x N„ 



(4) 



where A 7 ^- and are the number of /i and // + events, respectively, given by Eq. (|3|). Now 
N - is the number of yT events after taking care of the CID efficiency. All the quantities 
appearing in Eq. Q are functions of E^ and cos 6 p. 

Figure [3] shows the zenith angle distribution of events obtained before and after applying 
the reconstruction and CID efficiencies. Compared to Fig. [2j one can notice that the number 
of events fall sharply for the almost horizontal (cos 9 a ~ 0) bins. This can be understood by 



noting from Figure 1(c) that the reconstruction efficiency for muons falls as we go to more 
horizontal bins. The reason for this is the ICAL geometry where iron slabs and RPCs are 
stacked horizontally. Hence muons that are close to the horizontal direction cross fewer RPCs 
and mostly travel in the passive iron part of the detector. Identifying these tracks therefore 
becomes extremely difficult and the reconstruction efficiency of the detector falls. The CID 
efficiency also falls for the more horizontal bins and the result is that there are hardly any 
events for bins with —0.2 < cos 9^ < 0.2. 
Finally, the muon resolutions ge and a c 



5 are applied as follows: 



(5) 
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where (N®_)ij denotes the number of muon events in the i th E-bin and the j th cos#-bin after 
applying the energy and angle resolutions. Here E and cos 9 are the measured muon energy and 
zenith angle. The summation is over the true energy bin k and true zenith angle bin I, with E k 
and cos 9„ being the central values of the k th true muon energy and Z th true muon zenith angle 
bin. The quantities Kf and Mj are the integrals of the detector resolution functions over the 
bins of E and cos 9, the measured energy and direction of the muon, respectively. These are 
evaluated as: 

K '^=L iE ^7^ -4^ I • (6) 



and 



, , r^se Hj 1 / ( COS (9!, - cos6») 2 \ 
Mj(cos^) = / dcos9—= exp f~ 2 , (7) 



v27TO cos 0i i y 2(J COS0 ! 



where a E k and a cosd i are the energy and zenith angle resolutions, respectively, in these bins. We 
perform the integrations between the lower and upper boundaries of the measured energy (E^. 
and Ejji) and the measured zenith angle (cos 0£. and cos9h 3 )- For the extreme cos 9 bins, the 
bins are taken to be (— oo, -0.9) and [0.9, +oo) while integrating, and the events are assigned to 
the bins [-1, -0.9] and [0.9, 1], respectively. This ensures that no event is lost to the unphysical 
region and the total number of events does not change after applying the angular resolution. 
For E k < 1 GeV, the integrand in Eq. (6j) is replaced with the Landau distribution function, as 
the reconstructed energy distribution obtained from ICAL simulations |30| is specified in terms 
of this function. 
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Figure 4: Zenith angle distribution of \i 
after including energy and zenith angle resolution function. Here E and 9 are the measured 
energy and measured zenith angle, respectively. The error bars shown here are statistical. 
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Figure |4| shows the zenith angle distribution of fi~ events before and after folding in the res- 
olution functions. The angular dependence seems to get only slightly diluted. This is attributed 
to the good angular resolution of the detector. 

Table [3] shows the total number of muon events with measured energy range 0.8-10.8 GeV 
at various stages of the analysis for an exposure of 50 kt x 10 years. Note the sharp fall in 
statistics due to the reconstruction efficiencies. The reconstruction efficiencies are particularly 
poor for the near-horizontal bins where the reconstruction of the muon tracks is very hard. The 
small increase in the number of events after applying the energy resolution function is due to 
the spill-over of events from the low-energy part of the spectrum to measured energies greater 
than 0.8 GeV. The spillover to the energy bins with >10.8 GeV is comparatively small. The 
zenith angle resolution leaves the number of muon events nearly unchanged. 





A* 




Unoscillated 


14311 


5723 


Oscillated 


10531 


4188 


After Applying Reconstruction and CID Efficiencies 


4941 


2136 


After Applying (E, cos 9) Resolutions 


5270 


2278 



Table 3: Number of muon events produced in CC interactions at various stages of the analysis 
for an exposure of 50 kt x 10 years in the energy range 0.8-10.8 GeV. 



3.4 The x 2 Analysis 

For the x 2 analysis, we re-bin the data in ten energy bins between E = 0.8 — 10.8 GeV with 
bin width 1 GeV, and twenty cos#-bins between cos# = —1 to +1, with bin width 0.1. The 
simulated data for ICAL@INO is then statistically treated by defining the following x 2 '■ 



Xino 



(fi ) = min 



N E N cose 

EE 

i=l j=l 



pred 
ij 



n: 



obs\ 



2N°f s In 



^bs 



2 

k ' 



k=l 



with 



pred 



n: 



ij 



k=i j 



(8) 



(9) 



Here iV? red and N°f s are the expected and observed number, respectively, of \x events (N^_) 
in a given (E, cos9) bin, while Ne (=10) is the number of energy bins and N cos g (=20) is 
the number of zenith angle bins. N°^ s is calculated for a set of assumed "true value" of the 
oscillation parameters, listed in Table [2J Nf- is the predicted number of events for a given set 
of oscillation parameters without the systematic errors included. The systematic uncertainties 
are included via the "pull" variables one each for every systematic uncertainty a^. Here nfj 
is the change in the number of events in the (ij) th bin caused by varying the value of fc th pull 
variable by a^. For determining ir^, we have used a procedure similar to the one described 



in [39]. 

In this analysis we have considered the following five systematic uncertainties. We take 20% 
error on the flux normalization, 10% error on cross sections, and an overall 5% error on the total 
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number of events. In addition, we take a 5% uncertainty on the zenith angle dependence of the 
flux, and an energy dependent "tilt error" is included according to the following prescription. 
The event spectrum is calculated with the predicted atmospheric neutrino fluxes and then with 
the flux spectrum shifted according to 

**(£) = *o(£) (§-] ^ME) (l + 6\n^-) , (10) 

where Eq = 2 GeV and 5 is the la systematic tilt error, taken to be 5%. The difference between 
&s(E) and &o(E) is then included as the error on the flux. 

For each set of oscillation parameters, we calculate the x 2 separately for the \i~ and fi + data 
samples, and add them to obtain the total x 2 as 

xL = xLOO + xL(/^ + ) • (11) 

Since in this analysis we are mainly interested in constraining 623 and |Am 2 2 |, and the 
variation of 6*12 or lAm^l within the current error bars is observed not to affect the results, we 
take the value of these two parameters to be fixed to those given in Table [2] On the parameter 
sin 2 2#i3, we impose a prior to allow for the uncertainty in its current measurement : 

2 2 , f sin 2 2fli 3 (true) -sin 2 2fli3 V , 10 , 

X = Xino + , (12) 

V °sin 2 20 13 / 

where a sin 2 2 Q 13 is the current la error on sin 2 26*13, and is taken as 0.013 in our analysis. Of 
course, during the operation of INO, the error <r sin 2 2 g 13 will decrease, and within a few years, 
sin 2 26*13 may be considered to be a fixed parameter. 

4 Precision Measurement of sin 2 #23 and lAmf^l 

We start by presenting the reach of the IC AL for the parameters sin 2 6*23 and | Am| 2 1 separately. 
The true values of all parameters are given in Table [2J Note that we use the parameter Am| 2 ( 
Amfi - Am|i) instead of |Am|i|, the current limits on which are given in fl3~|. The x 2 values 



as functions of sin 2 6*23 and | Am 2 2 | are shown in figs. 5(a) and 5(b) , respectively. Note that the 
minimum value of x 2 vanishes, since MC fluctuations in the observed data have been reduced 
due to the scaling from an exposure of 50 kt x 1000 years. 
The precision on these parameters may be quantified by 

. . Pmax Pmin /-i n\ 

precision = , (Id) 

Pmax ~i~ Pmin 

where p m ax and p m in are the largest and smallest value of the concerned oscillation parameters 
determined at the given C.L. from the atmospheric neutrino measurements at ICAL for a given 
exposure. We find that after 5 years of running of this experiment, ICAL would be able to 
measure sin 2 6*23 to a precision of 20% and |Am| 2 | to 7.4% at la. With 10 years exposure, 
these numbers improve to 17% and 5.1% for sin 2 #23 and |Am 2 2 |, respectively. The precision 
on sin 2 6*23 is mainly governed by the muon reconstruction efficiency and is expected to improve 
with it. It will also improve as the systematic errors are reduced. If the flux normalization 
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Figure 5: The panel (a) shows the x 2 as a function of sin 2 623 for lAmf^l = 2.242 x 10 3 eV 2 
and sin 2 #23(true) = 0.5 . The panel (b) shows the x 2 as a function of jAm^l for sin 2 623 = 0.5 
and I Am|,| (true) = 2.242 x 10" 3 eV 2 



error were to come down from 20% to 10%, the precision on sin 2 623 would improve to 14% for 
10 years of exposure. Reducing the zenith angle error from 5% to 1% would also improve this 
precision to ~ 14%. On the other hand, the precision on |Am| 2 | is governed by the ability of 
the detector to determine the value of LjE for individual events accurately. This depends on 
the energy- and cos 9- resolution of the detector. 

A few more detailed observations may be made from the x 2 plots in fig. [Hj From fig. 5(a) one 
can notice that the precision on 623 when it is in the first octant (sin 2 823 < 0.5) is slightly better 
than when it is in the second octant (sin 2 623 > 0.5), even though the muon neutrino survival 
probability depends on sin 2 2023 a t the leading order. This asymmetry about sin 2 623 = 0.5 stems 
mainly from the full three-flavor analysis that we have performed in this study. In particular, we 
have checked that the non-zero value of Amfi is responsible for the asymmetry observed in this 



figure. On the other hand, x asymmetry about the true value of [Am^l observed in Fig. 5(b) 
is an effect that is present even with a two-flavor analysis. 

The precisions obtainable at the IC AL for sin 2 #23 and | Am| 2 1 are expected to be corre- 
lated. We therefore present the correlated reach of ICAL for these parameters in figs. |6(a)| and 
[6(b)] These are the main results of our analysis. As noted above, our three-neutrino analysis 
should be sensitive to the octant of #23- Therefore we choose to present our results in terms of 
sin 2 623 instead of sin 2 26*23. Though the constant-x 2 contours still look rather symmetric about 
sin 2 6*23 = 0.5, that is mainly due to the true value of sin 2 6*23 being taken to be 0.5. The values 
of sin 2 #23 away from 0.5 would make the contours asymmetric and would give rise to some 
sensitivity to the octant of 823, as we shall see later. 
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Figure 6: The precision reach expected at ICAL in the sin 2 823 — |Am 2 2 | plane at various 
confidence levels. The black(broken), blue(dotted) and red(solid) lines show 68%, 90% and 99% 
C.L contours. The true values of sin 2 823 and |Am| 2 | used for generating data are shown by 
the black dots. The true values of other parameters used are given in Table [2} Panel (a) is for 
five-year running of the 50 kt detector while (b) is for ten years exposure. 



Note that in our analysis, we have taken cos# bins of width 0.1. The muon angle resolution 
in ICAL is however better than 1° for almost all values of the zenith angle. The reason we had 
to take such wide angle bins is because of limited statistics. In order to study the impact of 
our choice of binning on the precision measurements of sin 2 823 and |Am| 2 | in ICAL, we reduce 
the cos 8 bin size to 0.05 and the size of the E bins to 0.5 GeV. In Fig. [7j we show the effect of 
taking these finer bins in E and cos 8 on the precision reach of the two parameters. The figure 
corresponds to 10 years of ICAL exposure. We see that the finer bins bring only a marginal 
improvement in the precision measurement of both sin 2 823 and | Am| 2 1 . For further results in 
this paper, we continue to use energy bins of width 1 GeV and zenith angle bins of width 0.1. 
The optimization with variable bin widths will be a part of future work. 

In fig. [8] we show the comparison of the precision reach on the atmospheric neutrino oscillatino 
parameters at ICAL with that obtained from other experiments currently. Note that here we use 
the parameter sm 2 26?23 instead of sirP'823 in order to enable a direct comparison. The blue and 
red lines show the expected sensitivity from atmospheric neutrino measurements at ICAL after 
5 years and the 10 years exposure, respectively. The green line is the 90% C.L. -allowed contour 
obtained by the zenith angle analysis of SK atmospheric neutrino measurements (shown at the 
Neutrino 2012 conference), while the pink line is the contour obtained by their L/E analysis [40] . 
The black line shows the 90% C.L. allowed region given by the combined analysis of the full 
MINOS data including 10.71 x 10 20 POT for the z^-beam, 3.36 x 10 20 POT for the z^-beam, as 
well as the atmospheric neutrino data corresponding to an exposure of 37.9 kt-years [4l]. 
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Figure 7: Effect of binning on the precision reach expected from the analysis of 50 kt x 10 years 
data of atmospheric neutrinos in ICAL. The black, blue and red lines show 68%, 90% and 99% 
CL contours. The solid lines show the contours for (E, cos#) bin widths of (1 GeV, 0.1) while 
the broken lines show the contours for (E, cos#) bin widths of (0.5 GeV, 0.05). 
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Figure 8: Comparison of the reach of ICAL@INO with the current results from other neutrino 
experiments. The black dot in the figure denotes the point where the ICAL data was generated. 
The true values of the other oscillation parameters are given in Table [2] 
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From fig. [8] it may be observed that, with 5 years of exposure, ICAL will be able to almost 
match the precision on |Am 2 2 | obtained from the SK L/E analysis currently. With 10 years 
data this will improve, though it will still not be comparable to the precision we already have 
from the MINOS experiment. Since the direction of neutrinos in MINOS is known accurately, 
their L/E is known to a greater precision and their measurement of |Am| 2 | is consequently more 
accurate. The precision of ICAL on sin 2 2023 i n 10 years may be expected to be comparable 
to what we currently have from SK. (Of course by the time ICAL completes 10 years, SK 
would have gathered even more data.) This precision is controlled to a large extent by the 
total number of events. We can see that the sensitivity of ICAL to sin 2 2^23 and |Am 2 2 | is 
not expected to surpass the precision we already have from the current set of experiments. In 
fact, the precision on these parameters are expected to improve significantly with the expected 
data from T2K [9,42 and NOi^A 21 , and ICAL will not be competing with them as far as 
these precision measurements are concerned. The ICAL data will however give complementary 
information on these parameters, which will significantly contribute to the improvement of the 
precision on the global fit. 

Earth matter effects in atmospheric neutrinos can be used to distinguish maximal from non- 



maximal #23 mixing and can lead to the determination of the correct #23 octant 43-45 . We 



show in Fig. [9] the potential of 10 years of ICAL run for distinguishing a non-maximal value 
of #23 from maximal mixing in the case where sin 26*23 = 0.90 (sin 2 <9 23 = 0.342, 0.658) and 
sin 2 2(9 23 = 0.95 (sin 2 6 23 = 0.388, 0.612). Note that the current 3cr allowed range of sin 2 2# 2 3 is 
(0.91, 1.0). The figure shows that, if the value of #23 is near the current 3a bound and in the first 
octant, then it may be possible to exclude maximal mixing to 99% C.L. with this 2-parameter 
analysis. If #23 is in the second octant, or if sin 2 2# 2 3 is larger than 0.9, the exclusion of the 
maximal mixing becomes a much harder task. 

Fig. [9] can also be used to quantify the reach of ICAL for determining the correct octant of 
#23, if the value of sin 2 2#23 is known. This can be seen by comparing the \ 2 value corresponding 
to the true value of sin 2 #23, but in the wrong octant, with that corresponding to the true value 
of sin 2 #23- We find that, for sin 2 2#23 = 0.9, i.e. just at the allowed 3a bound, the octant can 
be identified at >95% C.L. with 10 years of ICAL run if #23 is in the first octant. However if #23 
is in the second octant, the identification of the octant would be much harder: #23 in the wrong 
octant can be disfavored only to about 85% C.L.. The situation is more pessimistic if sin 2 2#23 
is closer to unity. 

The precision on |Am 2 2 | will keep improving with ongoing and future long baseline exper- 
iments. The inclusion of the information may improve the chance of ICAL@INO being able 
to identify deviation of #23 from maximal mixing and its octant to some extent, however that 
exercise is beyond the scope of this paper. 
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Figure 9: The projected reach in the sin 2 623 — |Am 2 2 | plane for four different non-maximal 
choices of 623- The black(broken), blue(dotted) and red(solid) lines show 68%, 90% and 99% 
C.L. contours for 10 years of 50 kt ICAL run. Note that we use normal hierarchy, and assume 
that it is already known. 



5 Summary 

In the analysis, we have calculated the projected reach of the ICAL experiment at INO for precise 
determination of the atmospheric neutrino parameters. We have used NUANCE simulated data 
on muon energies and directions using the resolutions and efficiencies obtained from a complete 
detector simulation by the INO collaboration. We incorporate the oscillation physics using the 
re- weighting method and calculate the values of x 2 m the sin 2 ^23 — I ^■ 7T? '32 1 pl ane using the pull 
method that takes into account systematic errors on (i) flux normalization (ii) cross sections 
(iii) neutrino flux tilt (iv) zenith angle dependence of flux and (v) overall systematic error. 

We present uncorrelated as well as correlated constraints on the value of sin 2 623 and | Am 2 2 1 
expected to be obtained after 10 years of running of the 50 kt ICAL. We find that the values 
of sin 2 623 and |Am 2 2 | may be determined at an accuracy of 17% and 5.1% respectively. The 
sensitivities with the data at ICAL only are not expected to be better than what we already 
have, indeed some of the other experiments in the next decade may do much better. However the 
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measurement at ICAL will be complementary and may be expected to contribute significantly 
towards the precision of parameters in a global fit. 

As far as the sensitivity to the 6*23 octant and its deviation from maximality is concerned, 
we find that 10 years of ICAL can exclude maximal mixing or the 623 value in the other octant 
to >95% C.L. only if the actual #23 is in the first octant and close to the current 3a lower 
bound to 99 % C.L.. Indeed, the octant identification seem to be beyond the reach of any single 
experiment in the next decade. 

This paper presents the first study on the reach of ICAL@INO for precision of atmospheric 
neutrino parameters, using the complete detector simulation. Note that in this analysis we have 
used only the information on muon events. However, the ICAL atmospheric neutrino experiment 
will also record and measure the hadrons associated with the charged current interaction of 
the muon- type neutrinos. Inclusion of this data-set into the analysis is expected to provide 
energy and angle reconstruction of the neutrino. This could lead to an improved sensitivity of 
the detector to the oscillation parameters. The analysis including the hadrons along with the 
muons is a part of the ongoing effort of the INO-ICAL collaboration. In addition, the updates 
and improvements in the muon momentum reconstruction algorithm as well as the optimization 
of our analysis procedure are likely to improve the results presented in this paper. 
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